Search for: All records

In this study, we apply optimal control theory to an immuno-epidemiological model of HIV and opioid epidemics. For the multi-scale model, we used four controls: treating the opioid use, reducing HIV risk behaviour among opioid users, entry inhibiting antiviral therapy, and antiviral therapy which blocks the viral production. Two population-level controls are combined with two within-host-level controls. We prove the existence and uniqueness of an optimal control quadruple. Comparing the two population-level controls, we find that reducing the HIV risk of opioid users has a stronger impact on the population who is both HIV-infected and opioid-dependent than treating the opioid disorder. The within-host-level antiviral treatment has an effect not only on the co-affected population but also on the HIV-only infected population. Our findings suggest that the most effective strategy for managing the HIV and opioid epidemics is combining all controls at both within-host and between-host scales. 

This study examines the interactions between healthy target cells, infected target cells, virus particles, and immune cells within an HIV model. The model exhibits two equilibrium points: an infection-free equilibrium and an infection equilibrium. Stability analysis shows that the infection-free equilibrium is locally asymptotically stable when R0<1. Further, it is unstable when R0>1. The infection equilibrium is locally asymptotically stable when R0>1. The structural and practical identifiabilities of the within-host model for HIV infection dynamics were investigated using differential algebra techniques and Monte Carlo simulations. The HIV model was structurally identifiable by observing the total uninfected and infected target cells, immune cells, and viral load. Monte Carlo simulations assessed the practical identifiability of parameters. The production rate of target cells (λ), the death rate of healthy target cells (d), the death rate of infected target cells (δ), and the viral production rate by infected cells (π) were practically identifiable. The rate of infection of target cells by the virus (β), the death rate of infected cells by immune cells (Ψ), and antigen-driven proliferation rate of immune cells (b) were not practically identifiable. Practical identifiability was constrained by the noise and sparsity of the data. Analysis shows that increasing the frequency of data collection can significantly improve the identifiability of all parameters. This highlights the importance of optimal data sampling in HIV clinical studies, as it determines the best time points, frequency, and the number of sample points required to accurately capture the dynamics of the HIV infection within a host. 

We used a previously introduced HIV within-host model with sensitive and resistant strains and validated it with two data sets. The first data set is from a clinical study that investigated multi-drug treatments and measured the total CD4+ cell count and viral load. All nine patients in this data set experienced virologic failure. The second data set includes a unique patient who was treated with a unique drug and for whom both the sensitive and resistant strains were measured as well as the CD4+ cells. We studied the structural identifiability of the model with respect to each data set. With respect to the first data set, the model was structurally identifiable when the viral production rate of the sensitive strain was fixed and distinct from the viral production rate of the resistant strain. With respect to the second data set, the model was always structurally identifiable. We fit the model to the first data set using nonlinear mixed effect modeling in Monolix and estimated the population-level parameters. We inferred that the average time to emergence of a resistant strain is 844 days after treatment starts. We fit the model to the second data set and found out that the all the parameters except the mutation rate were practically identifiable. 

In this paper we develop a four compartment within-host model of nutrition and HIV. We show that the model has two equilibria: an infection-free equilibrium and infection equilibrium. The infection free equilibrium is locally asymptotically stable when the basic reproduction number $$ \mathcal{R}_0 < 1 $$, and unstable when $$ \mathcal{R}_0 > 1 $$. The infection equilibrium is locally asymptotically stable if $$ \mathcal{R}_0 > 1 $$ and an additional condition holds. We show that the within-host model of HIV and nutrition is structured to reveal its parameters from the observations of viral load, CD4 cell count and total protein data. We then estimate the model parameters for these 3 data sets. We have also studied the practical identifiability of the model parameters by performing Monte Carlo simulations, and found that the rate of clearance of the virus by immunoglobulins is practically unidentifiable, and that the rest of the model parameters are only weakly identifiable given the experimental data. Furthermore, we have studied how the data frequency impacts the practical identifiability of model parameters. 

Understanding the epidemiology of emerging pathogens, such as Usutu virus (USUV) infections, requires systems investigation at each scale involved in the host–virus transmission cycle, from individual bird infections, to bird-to-vector transmissions, and to USUV incidence in bird and vector populations. For new pathogens field data are sparse, and predictions can be aided by the use of laboratory-type inoculation and transmission experiments combined with dynamical mathematical modelling. In this study, we investigated the dynamics of two strains of USUV by constructing mathematical models for the within-host scale, bird-to-vector transmission scale and vector-borne epidemiological scale. We used individual within-host infectious virus data and per cent mosquito infection data to predict USUV incidence in birds and mosquitoes. We addressed the dependence of predictions on model structure, data uncertainty and experimental design. We found that uncertainty in predictions at one scale change predicted results at another scale. We proposedin silicoexperiments that showed that sampling every 12 hours ensures practical identifiability of the within-host scale model. At the same time, we showed that practical identifiability of the transmission scale functions can only be improved under unrealistically high sampling regimes. Instead, we proposed optimal experimental designs and suggested the types of experiments that can ensure identifiability at the transmission scale and, hence, induce robustness in predictions at the epidemiological scale. 

Abstract Determining accurate estimates for the characteristics of the severe acute respiratory syndrome coronavirus 2 in the upper and lower respiratory tracts, by fitting mathematical models to data, is made difficult by the lack of measurements early in the infection. To determine the sensitivity of the parameter estimates to the noise in the data, we developed a novel two-patch within-host mathematical model that considered the infection of both respiratory tracts and assumed that the viral load in the lower respiratory tract decays in a density dependent manner and investigated its ability to match population level data. We proposed several approaches that can improve practical identifiability of parameters, including an optimal experimental approach, and found that availability of viral data early in the infection is of essence for improving the accuracy of the estimates. Our findings can be useful for designing interventions. 

In this paper, we present a multi-scale co-affection model of HIV infection and opioid addiction. The population scale epidemiological model is linked to the within-host model which describes the HIV and opioid dynamics in a co-affected individual. CD4 cells and viral load data obtained from morphine addicted SIV-infected monkeys are used to validate the within-host model. AIDS diagnoses, HIV death and opioid mortality data are used to fit the between-host model. When the rates of viral clearance and morphine uptake are fixed, the within-host model is structurally identifiable. If in addition the morphine saturation and clearance rates are also fixed the model becomes practical identifiable. Analytical results of the multi-scale model suggest that in addition to the disease-addiction-free equilibrium, there is a unique HIV-only and opioid-only equilibrium. Each of the boundary equilibria is stable if the invasion number of the other epidemic is below one. Elasticity analysis suggests that the most sensitive number is the invasion number of opioid epidemic with respect to the parameter of enhancement of HIV infection of opioid-affected individual. We conclude that the most effective control strategy is to prevent opioid addicted individuals from getting HIV, and to treat the opioid addiction directly and independently from HIV. 

In this paper, we introduce a novel multi-scale network model of two epidemics: HIV infection and opioid addiction. The HIV infection dynamics is modeled on a complex network. We determine the basic reproduction number of HIV infection, $$ \mathcal{R}_{v} $$, and the basic reproduction number of opioid addiction, $$ \mathcal{R}_{u} $$. We show that the model has a unique disease-free equilibrium which is locally asymptotically stable when both $$ \mathcal{R}_{u} $$ and $$ \mathcal{R}_{v} $$ are less than one. If $$ \mathcal{R}_{u} > 1 $$ or $$ \mathcal{R}_{v} > 1 $$, then the disease-free equilibrium is unstable and there exists a unique semi-trivial equilibrium corresponding to each disease. The unique opioid only equilibrium exist when the basic reproduction number of opioid addiction is greater than one and it is locally asymptotically stable when the invasion number of HIV infection, $$ \mathcal{R}^{1}_{v_i} $$ is less than one. Similarly, the unique HIV only equilibrium exist when the basic reproduction number of HIV is greater than one and it is locally asymptotically stable when the invasion number of opioid addiction, $$ \mathcal{R}^{2}_{u_i} $$ is less than one. Existence and stability of co-existence equilibria remains an open problem. We performed numerical simulations to better understand the impact of three epidemiologically important parameters that are at the intersection of two epidemics: $$ q_v $$ the likelihood of an opioid user being infected with HIV, $$ q_u $$ the likelihood of an HIV-infected individual becoming addicted to opioids, and $$ \delta $$ recovery from opioid addiction. Simulations suggest that as the recovery from opioid use increases, the prevalence of co-affected individuals, those who are addicted to opioids and are infected with HIV, increase significantly. We demonstrate that the dependence of the co-affected population on $$ q_u $$ and $$ q_v $$ are not monotone.